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SUMMARY 


In this paper it is demonstrated that COSMIC/ NAS TRAN may be used to analyze 
plate and shell structures made of symmetric composite laminates. Although general 
composite laminates cannot be analyzed using NASTRAN, the theoretical development 
presented herein indicates that the integrated constitutive laws of a symmetric 
composite laminate resemble those of a homogeneous anisotropic plate, which can be 
analyzed using NASTRAN. A detailed analysis procedure is presented, as well as an 
illustrative example. 


INTRODUCTION 


Composite laminates have been used in many engineering structures recently, 
and researchers in this field have developed many finite element programs [1]. To 
date, the composite laminate theory has not been incorporated into any NASTRAN 
plate element formulation [2]. 

The governing equilibrium equations for symmetric laminates resemble those for 
anisotropic plates [3]. Since NASTRAN possesses the capability to analyze aniso- 
tropic plates, symmetric laminates can be analyzed using NASTRAN if the equivalent 
anisotropic plate material properties can be obtained. 

The purpose of this study is to obtain equivalent anisotropic plate material 
properties for symmetric laminates consistent with the governing plate theory, and 
thus extend the capability of NASTRAN to include analysis of symmetric laminates. 
This theoretical formulation is given in the next section. Subsequently, a step- 
by-step analysis procedure and an illustrative example problem are presented. 


THEORY 


The classical laminate theory [4] will be used. A Cartesian coordinate system 
is assumed wit' :he x-y plane located at the midplane of the laminate. The normal 
coordinate is and the total laminate thickness is h. The laminate consists of 
perfectly bonded laminae. If h is considered small compared to the x and y dimen- 
sions of the laminate, the Kirchoff hypothesis may be assumed to be applicable. 

The displacements can then be written as 
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where Uo, Vo, and Wo are the midplane displacements in the x, y, and z directions, 
respectively. The strains at distance z from the midplane can be written as 
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The constitutive law at any point in the laminate may be written in the form 
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Relations between the stress and moment resultants and the stretching and 
curvature strains are obtained by integrating Equation (4) through the thickness; 
specifically, 
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where 
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Presently, the general form Equation (5) is not used for any plate element in 
NASTRAN; however, if the laminate is symmetric, then and thus the inplane 

extension and out of plane bending effects decouple. Equivalent anisotropic 
constitutive laws can be written for each effect. For in plane extension, the 
equivalent stress-strain law is 
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and for bending it is 
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The coefficients in Equations (7) and (8) are denoted as the equivalent anisotropic 
material properties. Since NASTRAN can be used for the stress analysis of a homoge- 
neous anisotropic plate, it can thus be used for symmetric laminate stress analysis. 


THE ANALYSIS PROCEDURE 


incorporating the aforementioned theory into NASTRAN for symmetric laminate 
stress analysis is not a straightforward task. The analysis procedure may be out- 
lined as follows: 


(1) The theoretical development presented above is used to calculate the coef- 
ficients of Equations (7) and (8). The results are inputted to NASTRAN 

by using MAT 2 cards. 

(2) To analyze plate or shallow shell structures, plate elements CTRIA1 and/or 
CQUAD1 are used, with associated MAT2 cards. Output will include grid 
point displacements and the midplane stretching and curvature strains of 
each element. 


(3) 

(4) 


The strains in each layer of the element are calculated using Equation (2) 
and the stresses in the layer are calculated using Equation (4). 

The stresses in each layer are transformed to the principal material 
directions and the failure criterion [5] 
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is used for ply failure evaluation. If the value of F exceeds unity, 
that layer is failed. 


EXAMPLE PROBLEM 


A Kevlar 49*/Epoxy cylindrical shell with [-60/0/60] 3S layup is used to illus- 
trate the analysis procedure. The geometry and finite element mesh are shown in 
Figure 1. The mechanical properties in principal material directions 1 and 2 of a 
unidirectionally reinforced lamina are 
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E 1 1 

= 7.6 X 

10 10 Pa 

E 2 2 

= 5.5 X 

10 9 Pa 

Gl2 

= 2.1 X 

10 9 Pa 

V 12 

= 0.34 



where the 1 and 2 directions are parallel and normal to the fiber direction, 
respectively. The ultimate strengths of the lamina are 

(5 1 ) T = 1-38 X 10* Pa 
(S j > c * -2. 76 X 10 8 Pa 

(5 2 ) t = 2.96 X 10 7 Pa 

(S 2 ) c = -1.38 X 10® Pa 
S 12 = 6.00 X 10 7 Pa 

where T and C denote tension and compression, respectively. The shell was sub- 
jected to 6395 Pa (1 psi) uniform pressure and its periphery was pinned. 


The equivalent homogeneous anisotropic plate constitutive law for in plane 
calculation was calculated as 
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The input to NASTRAN was made according to Step (2) of che procedure given above. 

The output displacement contours in the x, y, and z directions are shown in 
Figures 2, 3, and 4, respectively. Note that these deformations are not symmetric 
with respect to the centerlines of the shell . This is due to the presence of the 
D 16 and D 26 terms. 

Table 1 summarizes the layer stresses in principal material directions for a 
typical element. The small values of F indicate that no ply failures occur. 


*Kevlar 49 aramid fibers, manufactured by E. I. duPont de Nemours & Co., Inc. 
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CONCLUSIONS AND SUGGESTIONS 


NASTRAN has been demonstrated to be a feasible tool for the stress analysis of 
symmetric composite laminates, though equivalent anisotropic material properties 
and layer stress calculations must be performed outside the NASTRAN framework. To 
increase the NASTRAN analysis capability to unsymmetric laminate and to simplify 
the analysis procedure, it would be worthwhile to incorporate the newly developed 
composite laminate finite element programs into N^iSTRAN in order to analyze the 
ever-increasing number of composite laminate structures. 
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TABLE 1. LAMINATE STRESSED N PRINCIPAL MATERIAL 
DIRECTIONS AND -LURE PARAMETER (F) OF 
ELEMENT 85 


T 12 (Pa) 



-2.048 X 10 s 
-8.963 X 10 s 
1.144 X 10 6 
-2.930 X 10 5 
-8.481 X 10 5 
1.186 X 10 6 
-3.820 X 10 5 
-7.998 X 10 s 
1.227 X 10 6 
1.234 X 10 6 
-7.515 X 10 5 
-5.288 X 10 s 
1.276 X 10 s 
-7.033 X 10 s 
-6.178 X 10 5 
1.317 X 10 5 
-6.536 X 10 5 
-7.033 X 10 s 









FIGURE 1 - SYMMETRIC LAMINATED COMPOSITE SHELL WITH ["-60/0/601 
LAY-UP AND 0 DEGREE IS ALONG THE X-AXIS L 
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CONTOUR 
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(ran) 

1 

-4.857 
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CONTOUR DISPLACEMENT 

NO. (mm) 
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-14.248 
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3.150 
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8.949 

11.849 


FIGURE 3 - Y -DISPLACEMENT CONTOUR PLOT 
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